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Abstract 

Despite more than a quarter century of research, the nature of the second-order phase transition in the heavy-fermion 
metal URu2Si2 remains enigmatic. The key question is which symmetry is being broken below this "hidden order" 
transition. We review the recent progress on this issue, particularly focusing on the thermodynamic evidence from 
very sensitive micro-cantilever magnetic torque measurements that the fourfold rotational symmetry of the underlying 
tetragonal crystal is broken. The angle dependence of the torque under in-plane field rotation exhibits the twofold 
oscillation term, which sets in just below the transition temperature. This observation restricts the symmetry of the 
hidden order parameter to the E + - or E -type, depending on whether the time reversal symmetry is preserved or not. 
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1. Introduction 

At Th = 17.5 K, URu2Si2 undergoes a second-order 
phase transition accompanied by the large anomalies 
in thermodynamic and transport properties 01 0, S 01 • 
Several remarkable features have been reported in the 
lower-temperature hidden order phase |5j. URu2Si2 ex- 
hibits the body-centered tetragonal crystal structure [see 
inset of Fig. [3, and to date no structural phase transi- 
tion is observed at 7V Tiny magnetic moment appears 
(Mo ~ 0.03//b) below Tj, [6], but it is by far too small 
to explain the large entropy released during the transi- 
tion and seems to have an extrinsic origin ||7l m ISD ■ The 
electronic excitation gap is formed at a large portion of 
the Fermi surface and most of the carriers (~ 90%) dis- 
appears ifToL [n, 12]. The gap is also formed in the 
incommensurate magnetic excitations below the hid- 
den order transition foil . With application of pres- 
sure, the large-moment antiferromagnetic state with the 
wave vector Q c = (0,0,1) = (1,0,0) appears, which 
is separated from the hidden-order phase by the first- 
order phase transition [see Fig.Q] [5, 9]. The quantum 
oscillation experiments [14] show that the Fermi sur- 
face of the hidden-order phase is quite similar to that 
of the pressure-induced antiferromagnetic phase. These 



results suggest that the hidden order and antiferromag- 
netism are nearly degenerate. 

In general, a second-order phase transition causes a 
change in various type of symmetries, such as crys- 
tal, rotational, gauge and time reversal symmetries. 
An order parameter is introduced to describe the low- 
temperature ordered phase with reduced symmetries. 
Therefore the key to the nature of hidden order lies in 
understanding which symmetry is being broken. Sev- 
eral microscopic models, including multipole order- 
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|, orbital currents 12311 and helicity order [24], 
have been proposed. However, in spite of the intensive 
experimental and theoretical studies, what is the gen- 
uine order parameter in the hidden order phase is still 
an open question. 



2. Searching for fourfold rotational symmetry 
breaking 

2.1. Information from the superconducting state em- 
bedded in the hidden-order phase 

Inside the hidden-order phase, the unconventional su- 
perconductivity appears below T c = 1.4 K. When the 
antiferromagnetic phase is induced by the pressure, the 
superconductivity disappears [see Fig. [TJ: i.e. the super- 
conductivity can coexist with the hidden order, but not 
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Figure 1: Schematic temperature-pressure phase diagram of URU2S12. 
At ambient pressure, the second-order phase transition from the para- 
magnetic, heavy Fermi liquid metallic state to the hidden order (HO) 
state occurs at T/, = 17.5 K. At lower temperature below T c = 1.4K 
superconducting (SC) state appears. Above the critical pressure 
(~ 0.75 GPa), the large-moment antiferromagnetic (AF) state sets 
in through the first-order phase transition. The inset illustrates the 
schematic crystal structure (large white spheres, middle red spheres, 
small blue spheres are U, Ru, and Si atoms, respectively) with the spin 
arrangements in the antiferromagnetic phase. 



dynamical susceptibility at Q c and hence the pairing in- 
teraction becomes large and repulsive at Q c , leading to 
the sign change of the superconducting order parame- 
ter A(k) = -A(k + Q c ) B28I1 . The disappearance of su- 
perconductivity in the pressure-induced antiferromag- 
netic phase implies that the Q c spin fluctuations are sup- 
pressed by the antiferromagnetic ordering. Recent the- 
oretical calculations also point to the chiral cZ-wave su- 
perconductivity produced by the antiferromagnetic fluc- 
tuations J29h . 

The + signs in Eq. ([]]) stems from the degenerate na- 
ture of the two states under the tetragonal symmetry 
in the system. The two states are characterized by the 
phase winding directions of the order parameter in the 
ab plane as shown in Fig. [2] Recent lower critical field 
H C \(T) measurements detect some unusual anomaly be- 
low the superconducting transition B30TI . From this ob- 
servation, a possibility of the phase transition inside the 
superconducting state due to the degeneracy of the chi- 
ral (f-wave has been pointed out. Namely, if the hidden 
order state breaks the fourfold rotational symmetry, the 
superconducting states with two different chiralities (the 
states with + and - sign in the Eq. ([TJ form) may have 
different transition temperatures. Although the origin 
of the H C \{T) anomaly should be scrutinized, this has 
triggered the search for the in-plane fourfold symmetry 
breaking in the hidden order phase. 



with the antiferromagnetic order. Therefore the infor- 
mation on the superconducting state is also important 
for understanding the hidden order. The thermal con- 



ductivity 11121 12511 and specific heat [26] measurements 
using ver y cl ean single crystals, which become available 
recently 112711 . consistently show that the field depen- 
dence of these quantities has strong anisotropy. From 
these results, both suggest the existence of point nodes 
along the c axis in the superconducting order parameter. 
The symmetry analysis for the spin-singlet states leads 
to the chiral d-wave symmetry with the form IU2I1 



k, 
sin — c 
2 
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-a + i sin — -a 



(1) 



that has the zero points along the F and X lines in the 
Brillouin zone as well as the horizontal line nodes as 
shown in Fig. [2] This symmetry of superconducting or- 
der parameter A(k) has sign change at the points con- 
necting with the antiferromagnetic vectors Q c = (0, 0, 1) 
and (1,0, 0), which are equivalent for the body-centered 
tetragonal structure. This sign-changing order parame- 
ter is consistent with the superconductivity mediated by 
the antiferromagnetic spin fluctuations. The commensu- 
late antiferromagnetic fluctuations lead to the enhanced 



2.2. Magnetic torque as a sensitive probe for the rota- 
tional symmetry breaking 



Among experimental probes, the magnetic torque t = 
HoMV x H is a particularly sensitive probe for detect- 
ing the magnetic anisotropy, where V is a sample vol- 
ume and M is the induced magnetization. The torque 
is also a thermodynamic quantity which is the angle 
derivative of free energy. The torque in tiny crystals can 
be measured by using the micro-tip cantilever J3I, 321. 
which has a sensitivity orders of magnitude larger than 
the commercial SQUID magnetometer. This technique 
has been widely used for the de Haas-van Alphen quan- 
tum oscillation measurements at high magnetic fields or 
for determination of the anisotropy parameter in lay- 
ered superconductors. We use this technique to mea- 
sure in-plane anisotropy of magnetic susceptibility in 
the hidden-order phase 13311 . The torque measurements 
in magnetic fields H = (Hcos(/>, //sin0, 0) rotating 
within the tetragonal ab plane in URu2Si2 (where <p is 
the angle from the [100] direction) provide a stringent 
test whether the hidden order parameter breaks the crys- 
tal fourfold symmetry. In such a geometry r{<p, T, H) 
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Figure 2: Schematic arrangements of the Brillouin zone for the body-centered tetragonal structure of URu2Si2 (left). The superconducting order 
parameter with chiral d-wave symmetry described by Eq. Q has horizontal nodes in the planes shown by the red arrows as well as point nodes 
along the X and T lines indicated by the blue arrows. A cross sectional cut at k, = n/c including two neighboring Brillouin zones is also depicted 
(right) with the phase winding directions (circular arrows) for one of the degenerate chiral rf-wave order parameters. The other order parameter 
has the opposite phase winding directions in each position. At the two points connected by the antiferromagnetic wave vector Q c (green arrows), 
which is (0, 0, 1) (left) or (1, 0, 0) (right), the superconducting order parameter changes sign. 



has a twofold oscillation with respect to 0-rotation 

= \^0 Hl V [(Xaa ~ Xbb) Sin 2 <P ~ 2 Xab COS 2(f)] , (2) 

where the susceptibility tensor Xij is given by M,- = 
2^ 7 // ; [seeFig.Ha)]. 

i 

It should be stressed that in a system holding the 
tetragonal C\ symmetry, should be zero because 
Xaa = Xbb and Xab = 0. Even when we consider non- 
linear (higher order) susceptibility, the C\ symmetry al- 
lows only fourfold oscillation with respect to </>-rotation, 
and cannot lead to the twofold component. Finite val- 
ues of T20 may appear if a new electronic or magnetic 
state emerges that breaks the tetragonal symmetry. In 
such a case, rotational symmetry breaking is revealed 
by Xaa + Xbb or Xab + depending on the direction of 
the orthorhombic anisotropy. 

It has been reported by several groups l34l |9[] that 
URu2Si2 single crystals may contain ferromagnetic 
(FM) impurities. Indeed, some of our crystals also show 
the sign of the FM impurities: In these crystals, the 
torque curve in rotating field H — (H sin 9, 0, H cos ff) 
within the ac plane shows the hysteresis due to the FM 
impurities when the angle 6 from the c axis is close to 
90°. We have carefully selected small pure crystals, in 
which no such hysteresis is observed within the experi- 
mental resolution of ~ 10~ 14 Nm(~ 0.01% of the ampli- 
tude \A^b\ of the oscillation due to the difference between 
Xcc and^ flfl ). 

We also carefully aligned the magnetic field within 
the plane by using a system with two superconducting 



magnets generating H in two mutually orthogonal di- 
rections and a 3 He cryostat equipped on a mechanical 
rotating stage with a minimum step of 1/500 deg ll35ll . 
By computer controlling two magnets and the rotating 
stage, we were able to rotate H continuously within the 
ab plane with a misalignment less than 0.02 deg from 
the ab plane. With this technique, we were able to ob- 
serve no twofold oscillations in t(</>) measurements in 
the paramagnetic phase above 7), up to 40 K. 

3. Thermodynamic evidence for fourfold symmetry 
breaking 

3.1. Twofold oscillation in the in-plane magnetic torque 

Figure|3lb) displays the evolution of the in-plane 
torque t(0 ) at l^o^l = 4 T on entering the hidden order 
phase [33]. At all temperatures, torque curves are per- 
fectly reversible with respect to the field rotation direc- 
tions. Moreover, the torque curves remain unchanged 
for field cooling conditions with different field angles. 
t(0) can be decomposed as t = T20 + T44 + + ••; 
where T2 n <p = A-in sin 2n(<p - <po) is a term with 2n-fold 
symmetry with n = 1,2, • ■ •. In the middle and lower 
panels, the twofold and fourfold components obtained 
from the Fourier analysis are displayed. 

What is remarkable is that the twofold oscillation 
is distinctly observed in the hidden order state, while 
it is absent in the paramagnetic phase at T = 18 K 
slightly above T/,, as shown in the middle panels of 
Fig. Ob). Indeed, the torque curves shown by the up- 
per panels apparently become asymmetric with respect 
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Figure 3: (a) Schematic configuration of the magnetic torque measurements for afo-plane field rotation by using the micro-cantilever technique. 
The magnetic field H (blue arrow) induces the magnetization M (green arrow) in the URu2Si2 crystal. The torque along the c axis (red arrow) can 
be detected by the change in the piezo resister, which is measured by the bridge configuration, (b) Upper panels show raw magnetic torque curves 
as a function of the azimuthal angle <f> at several temperatures. All data are measured at \hqH\ = 4 T. Middle and lower panels show twofold ri^, 
and fourfold components of the torque curves which are obtained from the Fourier analysis. 



to 90° rotations below 7/,. We note that the fourfold 
oscillations (and higher-order terms) observed for 
all temperatures arise primarily from the higher-order 
nonlinear susceptibilities [36]. We therefore focus our 
attention on the twofold oscillation. Figure^a) de- 
picts the temperature dependence of the amplitude of 
the twofold symmetry \A^\. The amplitude onsets pre- 
cisely at 7/, = 17.5 K [see inset of Fig. 0a)]. This result, 
together with no hysteresis in the torque curves, again 
rules out a possibility that very tiny ferromagnetic im- 
purity is an origin of twofold symmetry. Based on these 
results, we conclude that the amplitude of twofold oscil- 
lations is a manifestation of intrinsic in-plane anisotropy 
of the susceptibility. The T2,j term has the form cos 20, 
which indicates the emergence of non-zero Xab below 
7ft. This means 



*[H0] =Xaa + Xab **[H0] = Xaa ~ Xa. 



(3) 



as schematically shown in Fig.Hfb). 

The in-plane anisotropy that sets in precisely at 7/, 
indicates that the rotational symmetry is broken in the 
hidden-order phase. The temperature dependence of 
1^2*1 <* Xab near T h is fitted with Xab « c(T h - T) + 
c'(Th - T) 2 much better than with a purely quadratic fit 
[see inset of Fig.UJa)]. The leading 7-linear term near 
T/, is naturally expected from the standard Landau the- 
ory of second-order transition when Xab is described by 
the squared term of an order parameter 77 cc (T h - T) 1 ^ 2 . 
Thus, the observed in-plane anisotropy IXablXaa im- 



plies that the hidden-order parameter tj breaks the four- 
fold tetragonal symmetry. 

3.2. Domain formation due to the degenerate nature of 
nematicity 

The breaking of the C4 rotational symmetry in a 
tetragonal crystal structure is in general a hallmark of 
electronic nematic phases 13711. In such nematic phases 
that reduce C4 to C2 symmetry, there are two degenerate 
directions to which the orthorhombic anisotropy is char- 
acterized. In the present case, the [110] (or [110]) and 
[110] (or [110]) directions are energetically equivalent 
(corresponding to the positive and negative Xab), which 
naturally leads to the domain formation in the hidden 
order phase [see Fig. [3b)]. 

The domain formation results in the crystal size de- 
pendence of the twofold oscillation amplitude in the 
torque experiment. When the crystal size is much larger 
than the domain size, the volume difference of the two 
different domains can be negligibly small compared 
with the total volume. Then the contributions of the 
two domains to the twofold oscillations in the torque 
have out of phase, which leads to the cancellation of the 
oscillations. Indeed, the torque measurements in a big- 
ger crystals show that the normalized signal becomes 
smaller, and in the mm-sized crystals the bulk suscep- 
tibility measurements cannot detect any difference be- 
tween ^[110] and ^-[1 10], as shown in Fig.Ua). Our 
results can be understood if the domain size is of the 
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Temperature (K) 

Figure 4: (a) Temperature dependence of the amplitude IA2J of the twofold oscillations in r((p) divided by the sample volume V for two crystals #8 
(red circles) and #7 (blue circles) with dimensions indicated in the figure. The corresponding susceptibility difference \2x a b\ = l^[H0] -^[110]| 
normalized by Xaa is indicated in the right axis. The susceptibility difference measured by the SQUID magnetometer for the large crystal #0 is also 
plotted (gray circles), (b) Schematic U atom arrangement (gray circles) in the ab plane. The emergent nonzero Xab is indicative of non-equivalent 
susceptibilities along [110] and [110] directions (arrows). 



order of tens of microns, which may explain the diffi- 
culties in observing this effect by using the conventional 
techniques. The fact that the torque curves remain un- 
changed for field cooling conditions with different field 
angles implies that formation of such domains is pre- 
dominantly determined by and strongly pinned by the 
underlying crystal conditions, such as internal stress or 
disorder. 

The magnitude of the susceptibility difference 
^[1 10] -;f[110] = Ixab deduced from the twofold os- 
cillations in the smallest crystal used in the torque mea- 
surements is as large as 10% of the in-plane suscepti- 
bility [see Fig. Ufa)]. We stress that even in the small 
crystals such a large signal cannot be accounted for by 
some unknown effect of the surface layer, whose vol- 
ume is of the order of 10~ 4 of the total volume. It is 
therefore unambiguous that the signal comes from the 
bulk of URu 2 Si 2 . 

3.3. Field dependence of the twofold oscillations 

In Figs.|5ja) and (b) we depict the amplitude of 
the twofold oscillation in t(^) divided by the field 
A2<t,/{fioH) as a function of H. The same plot for 
A2e/(f*()H) obtained from the ac-plane rotation is also 
shown in Fig.f5J a )- While a conventional paramagnetic 
//-linear dependence is seen in A2e/(f*oH), / '(poH) 
shows a highly unusual magnetic response. We note that 



the magnitude of A 20 /(jjqH) is of the order of a few G, 
which considerably exceeds the upper limit of remanent 
difference between the in-plane and c axis magnetiza- 
tion. This is another indication that any possible rema- 
nent effects of FM impurities cannot be associated with 
the observed twofold symmetry in the in-plane <p scans. 
We note that this magnitude of A2<t,/(fJ.oH) corresponds 
to the difference of the magnetizations along [110] and 
[110] directions M[110] - M[110] and hence is a mea- 
sure of the distribution of in-plane magnetic induction. 
We infer that this magnitude (~ 1G at 4T) is con- 
sistent with the anomalous broadening of the nuclear- 
magnetic -resonance (NMR) spectrum observed in the 
hidden order phase for H ± c n la, 13811 . which is also 
determined by the field distribution. 

The result of non-trivial field dependence of 
A2$/(fJ.oH) indicates that the in-plane magnetizations 
M[110](//) and M[110](//) have different (non-linear) 
field dependencies. The non-vanishing A2$/(fJ.oH) as 
jUo// — > T also suggests that the newly found in-plane 
anisotropy is not induced by magnetic field, but is likely 
to be an intrinsic property of the (zero-field) ground 
state of the hidden-order phase. Such a non-vanishing 
behavior in the zero-field limit may be associated with 
possible time reversal symmetry breaking. Further mea- 
surements sensitive at low fields are important to ad- 
dress this issue. 
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Figure 5: (a) Magnetic field dependence of the amplitudes of the twofold oscillations A2e/(MoH) for the oc-plane rotation and Aj^ldiaH) for the 
a/)-plane rotation, (b) Expanded view for Ao^KiiqH). 



4. Constraints on the hidden order parameter 

The broken fourfold rotational symmetry imposes 
strong constraints on theoretical models which att emp t 
to explain the hidden order 1I3I |40i l4ll I4H H H 14511 . 
Based on the observation of the twofold term in the 
magnetic torque, Thalmeier and Takimoto [41] made 
the Landau theory analysis for the symmetry calcifica- 
tion of the order parameters allowed in URu2Si2- The 
tetragonal symmetry group D\h of URu2Si2 has three 
nontrivial irreducible representations {A%, B\ and B2) 
and one doubly degenerate E representation. Each rep- 
resentation has the even (+) and odd (-) parities with 
respect to time reversal. Among these symmetries, they 
found that only the £-type order parameters can pro- 
duce the twofold oscillations consistent with the torque 
experiments. For degenerate E representation there are 
two possible solutions for the hidden order parameter 
characterised by the vector 77 = (j] a , r\b). The observed 
non-zero Xab indicates r\ oc (+1, +1), which leads to two 
domain types responsible for the phase of the twofold 
oscillations through the sign of r\ a r\b- The symmetry 
of the hidden order parameter is therefore restricted to 
E + - or £~-type, depending on whether the time reversal 
symmetry is preserved or not. 

Recently Fujimoto proposed the spin nematic order 
as a candidate for the hidden order B40I1 . which breaks 
the fourfold symmetry in the spin space. This may be 
categorized in the E + symmetry in the above represen- 
tation, but involves momentum dependence of the order 
parameter that may have nodes. He found that this or- 
der parameter can explain the temperature dependence 
of the Xab obtained by the torque experiments in a semi- 



quantitative way. 

More recently, Ikeda et al. used ab initio calculations 
to compare the multipole fluctuations in URu2Si2 114311 . 
From the complete set of calculations for all the allowed 
symmetries, they conclude that the £T-type rank-5 do- 
triakontapole order is the most stable order that emerges 
at low temperatures in URu2Si2- They also found that 
the antiferromagnetic state is nearly degenerate in en- 
ergy, which can naturally explain the pressure induced 
phase transition from the hidden-order phase to the an- 
tiferromagnetic state (Fig.Q} as a spin-flop transition 
of pseudo-spins from in-plane to out-of-plane direction. 
The same order has also been proposed by Rau and Kee 
Il44ll . who found that the calculated temperature depen- 
dence Xab in this ordered state is consistent with the ex- 
periment. Another £~-type order parameter, which has 
some similarlity to the spin nematic order, has been pro- 
posed by Chandra et al. i45ll . who named as a 'hastatic' 
(Latin: spear) order. Such E~ orders break the C4 sym- 
metry as well as the time reversal symmetry. Thus the 
experimental verifications on the time reversal symme- 
try breaking are highly desired. 



5. Concluding remarks 

We have reviewed several recent advances on the ex- 
perimental and theoretical studies of URu2Si2, focusing 
on the magnetic torque experiments that indicate the 
fourfold rotational symmetry breaking in the hidden- 
order phase. Electronic states that break crystal ro- 
tational symmetry without showing ordered moments 
have also been suggested in various strongly correlated 
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electron systems 13711 . including Sr3Ru207 I46H . high- 
T c cuprates [47], and iron pnictides [48], which are dis- 
cussed in terms of stripe or nematic orders. The finding 
of a directional electronic state in a heavy-fermion ma- 
terial is another example of such exotic order in corre- 
lated matter. The broken C\ symmetry imposes strong 
constraints on the symmetry of the hidden order param- 
eter to the two-component E representation. The deter- 
mination of the order parameter requires experimental 
input as to whether time reversal symmetry is also bro- 
ken or not. Information on the electronic structure in 
the hidden order phase is another important factor in the 
identification of the hidden order parameter. Recent cy- 
clotron resonance experiments allow the full determina- 
tion of the effective masses on the main Fermi surface 
sheets [49]. Studies of the Fermi surface structure for 
the proposed order parameters are a promising way to 
resolving the long-standing hidden order mystery. 
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